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ABSTRACT 

The millimeter/submillimeter wavelength polarization of Sgr A* is known to be variable in both 
magnitude and position angle on time scales down to a few hours. The unstable polarization has 
prevented measurements made at different frequencies and different epochs from yielding convincing 
measurements of Faraday rotation in this source. Here we present observations made with the Sub- 
millimeter Array polarimeter at 227 and 343 GHz with sufficient sensitivity to determine the rotation 
measure at each band without comparing position angles measured at separate epochs. We find the 
10-epoch mean rotation measure to be (—5.6 ± 0.7) x 10^ rad m~^; the measurements are consistent 
with a constant value. We conservatively assign a 3a upper limit of 2 x 10^ rad m~^ to rotation 
measure changes, which limits accretion rate fluctuations to 25%. This rotation measure detection 
limits the accretion rate to less than 2 x lO^*" Mq yT~^ if the magnetic field is near equipartition, 
ordered, and largely radial, while a lower limit of 2 x 10~^ Mq yr~^ holds even for a sub-equipartition, 
disordered, or toroidal field. The mean intrinsic position angle is 167°±7° and we detect variations of 
31^9^ degrees. These variations must originate in the submillimeter photosphere, rather than arising 
from rotation measure changes. 
Subject headings: black hole physics 
interferometric 
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1. INTRODUCTION 

The linear polariz ation of Sgr A* was first detected by 
lAitken et al.l ()200Cll ) above 100 G Hz, after un successful 
searches at lower frequencies fe.g. lBower et al...l99 9a,b). 
Subsequent interferometer observatio ns have shown the 
polarization t o vary in position ang le (|Bower "eiral][2005h 
and fraction (jMarrone et al.l [20061 : hereafter M06), with 
variability occurring on timescales comparable to those of 
the previously observed total intensity variations (MOO). 
The variability may be intrinsic to the source or due to 
propagation effects, but the short timescales suggest pro- 
cesses very close to the black hole are responsible and 
thus polarization should be a useful tool for the study of 
Sgr A*. 

One aspect of the polarization that has yet to be ex- 
ploited is its variation with frequency, particularly the 
frequency-dependent orientation change known as Fara- 
day rotation. The presence of polarization was used im- 
mediately after its detection to argue that the infalling 
plasma must be tenuous and the mass accretion rate 
low (less than 10~^ Mq yr~^), because larger accre- 
tion rates would depolarize the emission through extreme 
Faraday rotation gradients (l Ouataert fc Gruzinov|[2000al : 
IXgol 2000). Constraints have also been deri ved through 
careful exa mination of the spectrum (e.g. iMelia et al.l 
[2000l[200lh . However, further progress on determining 
the accretion rate and its variability has awaited mea- 
surement of the Faraday rotation measure (RM), a dif- 
ficult task because of significant variability of Sgr A* 
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and the diminished effect of rotation at high frequencies. 
To date there have been three RM determination s from 
non-simultaneous obs ervations (|Bower et al.l 120031 : lM06l : 
iMacquart et al.ll2006t). but non e of them is robust. Most 
recentlv. iMacauart et al.l (120061 ) estimated the RM to be 
—4.4 X 10^ rad m~^ from their 83 GHz polarization data 
and all previous data. However, their analysis allows a 
much lower RM due to the 180° degeneracy of the polar- 
ization position angle (see § [3]) . 

Simultaneous measurements are a much more secure 
way to determine the RM because they are insensitive to 
source variability. Such data can be used to measure the 
RM, the intrinsic polarization direction, and the variabil- 
ity of each. Isolating these changes is crucial to under- 
standing the conditions in the inner accretion flow; the 
RM can be used to place upper and lower limits on the 
accretion rate, while variations in the polarization and 
RM with frequency and time can be used to examine 
the structure of the flow. The instantaneous frequency 
coverage of available instruments has not yet been ade- 
quate to show Faraday rotation in a single epoch, with 
the lowest upper limit (|RM| <7 x 10^ rad ni~^) from 
early Submillime ter A rray*^ (SMA) 345 GHz polarimet- 
ric observations ()M06f ). Here we report SMA measure- 
ments of the RM and intrinsic polarization changes of 
Sgr A*, the ffist observations to detect the RM with 
statistical significance and the only measurement made 
without resorting to comparisons between position an- 
gles measured at different epochs. The observations and 
calibration procedures are detailed in § [21 the RM and 
the variability of the RM and intrinsic polarization are 
examined in § [3l and we discuss the implications for the 
accretion rate and properties of Sgr A* in § HI 

"* The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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2. OBSERVATIONS 

The SMA made polarimetric observations of Sgr A* in 
2005 June and July. In the ten nights without weather 
or technical problems there are six (four) tracks in the 
230 GHz (345 GHz) band. The 230 GHz observations 
used a local oscillator (LO) frequency of 226.9 GHz, pro- 
viding simultaneous measurements of sidebands centered 
on 221.9 and 231.9 GHz, except for one observation made 
at LO=225.55 GHz. The higher frequency data use an 
LO frequency of 343 GHz. The weather was excellent 
in the usable nights, with a typical zenith opacity at 
225 GHz of 0.06 - 0.08, or 0.19 - 0.27 when scaled to 
343 GHz. All data were obtained with five to seven an- 
tennas in the compact configuration (7 — 70 meter base- 
lines), resulting in a synthesized beam width of 2'.'0x4'.'0 
at 230 GHz and l'.'6x3'.'2 at 345 GHz on Sgr A*. For 
the purposes of self-calibration and polarization extrac- 
tion, baselines shorter than 20fcA were removed to ex- 
clude the extended emission around Sgr A*. Gain cal- 
ibrators included J1743-038, J1733-130, J1924-292, and 
Ceres. Very similar polarization results were obtained 
when the calibrator gains were applied and when Sgr A* 
was self-calibrated. The flux density scale was derived 
from observations of planets and their moons, with an 
expected accuracy of better than 10%. Flux densities 
in each Stokes parameter were obtained through point- 
source fits to the visibilities. The observed properties of 
Sgr A* in all ten epochs are listed in Table [T] 

To obtain full p olarization informati on we used the 
SMA polarimeter ()M06l ; iMarrond [20061) . Measurement 
of polarization relies on precise determination of the 
fractional contamination ( "leakage" ) of each polarization 
state by the cross-handed polarization. Uncalibrated 
leakage contaminates the linearly polarized Stokes pa- 
rameters {Q and U) with /. Leakages were measured 
each night by observing polarized quasars (3C 279 or 
3C 454.3) over a large range of parallactic angle. The 
leakages are stable, approximately 1% in the upper side- 
band (USB) for each band, 4% for the 230 GHz lower 
sideband (LSB), and 2% for the 345 GHz LSB. The 
r.m.s. leakage variability is 0.3% at 230 GHz and 0.4% at 
345 GHz. There are three potential sources of variabil- 
ity: (1) real instrumental polarization changes between 
nights, (2) finite signal-to-noise ratio, and (3) polariza- 
tion leakage that is not constant across the sky, resulting 
in leakage determinations that depend on the hour angle 
coverage. The first of these may exist but is bounded 
by the small observed leakage variability. The second is 
expected to yield variations at the level of 0.1% in the 
345 GHz leakages and less at 230 GHz. We know that 
the third effect is also present, at a level of around 0.2%, 
due to the fact that cross-polarization introduced sky- 
ward of the Nasmyth relay mirror is rotated relative to 
the feed as a source is tracked in elevation, while the 
cross-polarization of the wave plate and subsequent op- 
tics is fixed. This small additional leakage imprints itself 
on the polarization of the calibrator and its effect on the 
leakages derived in calibration (under the assumption of 
stationary leakage) changes with the elevation coverage 
of the calibrator. Fortunately, the resulting polarization 
changes are small. Leakage errors of 0.3%— 0.4% con- 
tribute 0.2% fractional contamination to the instanta- 
neous polarization, significantly less when averaged over 



a full track due to the parallactic angle rotation. For 
Sgr A*, which we measure to be 5%— 10% polarized, this 
results in at most 1° of position angle error. We reduced 
the Sgr A* data using the average leakages at each fre- 
quency instead of the single-night values to confirm this 
estimate. More importantly for this purpose, the uncor- 
rected instrumental polarization is very nearly constant 
across the sidebands because properties such as the illu- 
mination pattern on the antenna change slowly with fre- 
quency. We find the imprint of this additional leakage on 
the calibrator to be very consistent between sidebands, 
so although the absolute position angle varies by up to 
1°, the inter-sideband difference varies only by 0.1°— 0.2°. 
We take this as our systematic error for position angle 
differences, although it is much lower than the thermal 
noise. Previous SMA 345 GHz observations (M06) do 
not include enough calibration data to place similar lim- 
its on the systematic errors, so we exclude them from the 
following analysis. 

3. ROTATION MEASURE AND INTRINSIC POLARIZATION 

Faraday rotation changes the observed polarization po- 
sition angle (x) as a function of frequency according to: 

X(^) = Xo + (1) 

where is the intrinsic position angle. The rotation 
measure (RM) is proportional to the integral of the elec- 
tron density and magnetic field component along the line 
of sight. From the observed LSB and USB position an- 
gles we derive a RM and xo for each observation. We 
plot X in both sidebands against for all ten epochs in 
Figure [1] along with the average RM and xo • The larger 
errors in Table [T] for 345 GHz band RMs and xos are due 
to the much smaller inter-sideband difference in at the 
higher frequency; the constraints on these quantities are 
therefore dominated by the 230 GHz data. The average 
RM from all 10 epochs is (-5.6 ± 0.7) x 10^ rad m'^, 
while the 230 GHz and 345 GHz points alone yield 
(-5.4±0.7)xl0'5 rad m'^ and (-13±5)xl05 radm'^, re- 
spectively, consistent within their errors. The ten single- 
night RM values are consistent with a constant value 
(Xr = 1.21, for 9 degrees of freedom). None of the 
230 GHz points deviates from the average by more than 
1.2 times its measurement error, while the largest de- 
viation at 340 GHz is 2.2a. We therefore place a con- 
servative upper limit of three times the averaged ti, or 
2 X 10^ rad m~^ on RM variability. This agrees with 
the small RM fluctuations inferred from the co nstancy 
of polarization at 83 GHz (jMacauart et al.ll2006f ). 

These observations also constrain the intrinsic polar- 
ization direction of Sgr A*. The average xo in our ten 
measurements is 167° ± 7°, or 162° ± 7° and 210° ± 21° 
from 230 and 345 GHz observations, respectively. The 
Xo values vary by more than our measurement errors pre- 
dict, suggesting intrinsic polarization changes. Assuming 
a constant xoj we obtain Xr = 2.8 (0.3% probability) us- 
ing all data points, or x^ = 3.9 (0.16%) for the 230 GHz 
points only. The scatter in the full data set suggests an 
intrinsic xo dispersion of 31^g* degrees, with very sim- 
ilar results obtained from just the 230 GHz points. It 
is interesting to note that the mean xo is ~80° different 
from the angle observed during IR flares (jEckart et al.l 
120061: iMever et al]|2006f) . A 90° position angle change 
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TABLE 1 

Polarization and Rotation Measure of Sagittarius A* 
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^ Statistical errors only. Overall flux density scale uncertainty is 10%. 



—4.4 X 10^ rad m^^ from the average x ^-t four frequen- 
cies over the last several years. However, their interpre- 
tation requires a 180° unwrapping of the 83 GHz po- 
sition angle and the non-wrapped position angle is not 
strongly excluded. They report a Xr of 2.1 for 3 de- 
grees of freedom (10% probability) for the non-wrapped 
fit (RM=— 1.9 X 10^ rad m~^). Furthermore, this xl re- 
lies on their poorly measured USB polarization (just two 
of five measurements detect polarization at 3a) and the 
standard error of the mean of just two measurements at 
216 GHz in order to provide additional degrees of free- 
dom. Discarding the former and using the 230 GHz vari- 
ability to estimate the 216 GHz variability increases the 
probability of the non-wrapped position angle to 26%. 

4. DISCUSSION 

The KM we observe is too large to be produced by 
material beyond the accretion radius of Sgr A* (approx- 
imately 1" or 0.04 pc). Using the density determined by 
iBaganoff et al.l (j2003l ). the RM in the dense inner 10" is 
just 8 X 10'^ rad m~^ assuming a 1 mG ambient field. 
RMs determined in nearby sources are as large as 70% 
of this estimate fe.g. lYuseip-Zadeh fc Morris! 1 19871 ). 

The RM can be used to determine the accretion rate 
(Af ) at small radii around Sgr A* if assumptions are 
made about the nat ure of the accretion flow. The pro- 
cedure is outlined in lM06l and assumes a power-law ra- 
dial density profile (n oc r~^) and an ordered, radial, 
equipartition-strength magnetic field. Our RM detection 
places upper limits that are 15% lower than those in their 
Figure 4. The important parameter for determining the 
upper limit is rjn, the radius at which electrons become 
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Fig. 1. — The position angles observed in each sideband on each 
of the ten epochs. The slope between the two sidebands at each 
epoch is proportional to the RM and the extrapolated intercept at 
the X axis is xo- The mean RM— xo fit is also plotted. Variability 
in xo is visible in the 230 GHz points (at right), which show slopes 
similar to the mean but are widely dispersed in x- 



has been predicted to occur near the spectral peak in 
theoretical models of the polarization, although the pre- 
cise fl ip frequency depen ds on the details of the models 
(|Agoll 120001: iMelia et al.ll2000l) . If the IR measurements 
trace the intrinsic polarization direction at short wave- 
lengths, rather than the possibly random polarization 
of the transient flare, this would be evidence for such a 
change above 345 GHz. We do not observe a flip between 
230 and 345 GHz. 

This measurement represents the first reliable deter- 
mination of the RM of Sgr A*, and the only measure- 
ment made from simultaneous observations at multiple 
frequencies an d therefore able to isola te source polariza- 
tion changes. iMacauart et al.l (|2006f ) derived a RM of 
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relativis tic and their RM contribution begins to be sup- 
pressed (|Quataert fc Gruzinovll2000a[ ). In general, shal- 
low density profiles approaching th at of the convection- 
dominated accretion flow (0 l /2: lNaravan et ani200Cll : 
iQuataert fc Gruzinovl l2000bf ) yi eld hotter centr al tem- 
peratures and thus larger rin (jQuataertI l2003f l. while 
steep profiles (/3 3/2), as in Bo ndi-like or advection- 
dominated accretion flows (ADAF; iNaravan fc YilfTQQSf ) 
are marginally relativistic to s mall radii and may have 
rin of a few rs . Simulation s (e.g. 'Hawley fc Balbusll2002l : 
llgumenshchev et al. 20dl [Pmet al. 2003 |) favor < 1. 
Those with published temperature profiles show rin be- 
tween 30 and 100 Schwarzschild radii (rs), yielding upper 
limits for M of 2 x 10""^ to 5 x IQ-^ Mp, yr' ^ 

Note that the assumptions of the |M06| form alism, 
which were taken from previous rel ated works (e.g. lMelial 
Il992t [Ouataert fc Gruzinovl l2000al) . are not constrained 
by existing observations. In particular, the assump- 
tion of equipartition-strength fields cannot be justified 
observationally, although simulations may show a ten- 
dency to reach equipartiti on fractions of a few percent 
(jlgumenshchev et aLll2003[ ). Magnetic fields that are a 
fraction e of the equipartition strength will raise the ac- 
cretion rate limits by e"^/^ (a factor of 10 for e — 3%). 
The assumption of a non-reversed field is also poten- 
tially suspect. The direction of position angle change 
with frequency, and thus the sign of th e RM, appears 
unchanged since the lAitken et al.l (|2000( l measurements 
in 1999. This is not natural for a turbulent accre- 
tion fiow, except along spe cial lines of sight (e.g., th e 
axial region in model A of llgumenshchev et al.l l2003f ) . 
iRuszkowski fc BegelmanI (|2002l ) suggested that the sta- 
bility of the sign o f circular polarization (CP) noted by 
I Bower et al.l (l2002f ) could be explained by Faraday con- 
version of linear to circular polarization in a highly re- 
versed field with a small directional bias. The bias field 
(Bb) naturally creates stability in the CP and RM signs 

if i?f, 3> Brms/^/N, where Brms is the random field and 
N the number of reversals. In this limit the integrated 
product of the electron density and parallel magnetic 
field, without regard to field reversals, would be large, as 
would the accretion rate, but the net RM would be small 
a nd originate in the bias field. The turbulent jet model 
of iBeckert fcFa lckc (20.02!) also uses a random field with 
a small bias to produce stability. If these models are cor- 
rect the accretion rate upper limits derived from our RM 
detection would no longer hold. 

Our RM also allows us to place lower limits on the ac- 
cretion rate; these are not subject to the above caveats 
since the uncertainties act to raise the minimum ac- 
cretion rate. If we take Hn to be around lOrg or 3rs 
(smaller r^n yields smaller lower limits, so this is conser- 
vative for hot fiows), we find that M must be greater 
than 1 - 2 X 10"^ Mq yr'^ or 2 - 4 x 10"^ Mq yr-\ 
respectively. If the field is toroidal, reversed, or sub- 
equipartition, these lower limits are raised and may pose 
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problems for very low- M models. Using the RM-Af^/^ 
scaling derived in lM06l . the observed limit on RM vari- 
ability (35%) limits accretion rate variability to 25% over 
two months. 

Our detection of significant variability in xo is the first 
clear separation of the effects of a variable RM and vari- 
able intrinsic polarization; disentangling these effects re- 
quires multi-frequency observations. In the optically thin 
limit Xo represents the intrinsic polarization direction of 
Sgr A* at all frequencies and x variations should fol- 
low Xo variations, with RM fluctuations potentially con- 
tributing a frequency-dependent component to the x dis- 
persion. However, comparing all previous measurements 
near 230 and 345 GHz (employing our RM to translate 
all observations to a single frequency), we find a disper- 
sion of 22° in 16 observ ations at 230 GHz and just 8° in 
11 epochs at 345 GHz (jAitken et all 120001: iBower et al] 
120031 120051 : lM06t this work). Although these two disper- 
sions scale by approximately consistent with RM 
fluctuations of approximately 2 x 10^ rad m^^, we con- 
sider this explanation unlikely. From equation ([1]), the 
proportionality between RM and x fluctuations should 
be if RM changes are responsible for x changes, with 
an increasing RM (less negative) expected for increasing 
X- In fact our 230 GHz data show the opposite trend, 
with the least negative RM values occurring on the days 
with the smallest x- The best-fit slope of the RM— X230 
relation is (— 3.4 ± 3.6) x 10^ rad m~^ deg"'^, different 
from the expected slope of +10 x 10"^ rad m^^ deg~^ 
by 4cr. Similar analysis is inconclusive in the 345 GHz 
data, which are less sensitive to RM. This confirms that 
we are seeing intrinsic polarization changes rather than 
RM fiuctuations. Moreover, the necessary RM changes 
would yield variability of 150° at 83 GHz , much larger 
than the 8° observed by iMacauart et all ()2006f ). The 
emission from Sgr A* appears to be optically thick be- 
low 300 — 400 GHz (Marrone et al., in prep) so xo may 
vary with frequency, contrary to the expectation at low 
optical depth. The observed spectrum of position angle 
variability between 83 and 345 GHz most likely repre- 
sents differences in the stability of the magnetic field ori- 
entation (manifested as xo) a-t the photospheres at these 
frequencies. If the Faraday screen is located at much 
larger radii than the submillimeter emission region, its 
variations are likely to be slower than those of the in- 
trinsic polarization; this allows the intrinsic changes to 
be isolated from putative RM fluctuations by examining 
short timescales, even without better knowledge of the 
RM stability. 



The authors thank the SMA staff, particularly Ken 
Young, for assistance with the polarimeter, and Avi Loeb 
and Ramesh Narayan for useful discussions. DPM ac- 
knowledges partial support from a Harvard University 
Merit Fellowship. 



Bower, G. C., Backer, D. C., Zhao, J.-H., Goss, M., & Falcke, H. 

1999a, ApJ, 521, 582 
Bower, G. C., Falcke, H., Sault, R. J., & Backer, D. C. 2002, ApJ, 

571, 843 

Bower, G. C., Falcke, H., Wright, M. C., & Backer, D. C. 2005, 
ApJ, 618, L29 



Faraday Rotation in Sagittarius A* 



5 



Bower, G. C, Wright, M. C. H., Backer, D. C, & Falcke, H. 1999b, 
ApJ, 527, 851 

Bower, G. C., Wright, M. C. H., Falcke, H., & Backer, D. C. 2003, 
ApJ, 588, 331 

Eckart, A., Schodel, R., Meyer, L., Trippe, S., Ott, T., & Genzel, 

R. 2006, A&A, 455, 1 
Hawley, J. F. & Balbus, S. A. 2002, ApJ, 573, 738 
Igumenshchev, I. V., Narayan, R., & Abramowicz, M. A. 2003, 

ApJ, 592, 1042 

Macquart, J.-R, Bower, G. C., Wright, M. C. H., Backer, D. C., & 

Falcke, H. 2006, ApJ, 646, Llll 
Maxrone, D. P. 2006, PhD thesis, Hajrvard University 
Marrone, D. P., Moran, J. M., Zhao, J.-H., & Rao, R. 2006, ApJ, 

640, 308 
Melia, F. 1992, ApJ, 387, L25 

Melia, F., Liu, S., & Coker, R. 2000, ApJ, 545, L117 



— . 2001, ApJ, 553, 146 

Meyer, L., Schodel, R., Eckart, A., Karas, V., Dovciak, M., & 

Duschl, W. J. 2006, A&A, 458, L25 
Narayan, R., Igumenshchev, I. V., & Abramowicz, M. A. 2000, 

ApJ, 539, 798 
Narayan, R. & Yi, I. 1995, ApJ, 452, 710 
Pen, U.-L., Matzner, C. D., & Wong, S. 2003, ApJ, 596, L207 
Quataert, E. 2003, Astronomische Nachrichten Supplement, 324, 

435 

Quataert, E. & Gruzinov, A. 2000a, ApJ, 545, 842 
— . 2000b, ApJ, 539, 809 

Ruszkowski, M. & Begelman, M. C. 2002, ApJ, 573, 485 
Yusef-Zadeh, F. & Morris, M. 1987, ApJ, 322, 721 



